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Abstract 

We investigate ferromagnetism in the periodic Anderson model with diagonal disorder. Using dynamical mean- 
field theory in combination with the modified perturbation theory, the disorder can be included in the calculation 
consistently, which turns out to be equivalent to the CPA method. Disorder generally reduces the Curie temper- 
ature and can for certain configurations completely suppress ferromagnetic order. This can be ascribed to the 
enhanced quasiparticle damping and the special structure of the density of states. PACS 
71.10.Fd, 71.28+d, 75.30.Md 
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The investigation of models for strongly cor- 
related electron systems has made significant 
progress in recent years due to the introduction 
of the dynamical mean-field theory ||, ^ which 
is based on the non-trivial limit of infinite spa- 
cial dimensions In this limit, the lattice 

self-energy becomes wave-vector independent and 
the problem can be mapped onto a single-site 
problem Q. Many questions concerning strongly 
correlated electron systems such as the Mott- 
Hubbard metal-insulator transition could be an- 
swered by this approach ^. Also important 
insight into the physics of band-ferromagnetism 
could be gained ||, |j. 

The treatment of disorder also simplifies in the 
limit of infinite spatial dimensions [ pO| : The well- 
known CPA method which has to be seen 
as the best single-site approximation for solving 
disorder problems [|l2j, becomes exact here. 

As first shown by Ulmke et al. , the dynam- 
ical mean-field theory therefore allows to investi- 
gate the interplay of disorder and strong electron 
correlations by taking into account both prob- 
lems on the same level of approximation (see also 

§, 0). 

In this paper we want to focus on the influence 
of disorder on band-ferromagnetism, in particular 
the ferromagnetic phase of the periodic Anderson 
model in the intermediate-valence regime [TtI 
[l9|. The periodic Anderson model is defined by 



its Hamiltonian: 



,(f)^(f) 



(1) 



Here, s^^ (fia) and st^ ifL) ^he creation 
and annihilation operators for a conduction elec- 
tron with Bloch vector k and spin a (a localized 



electron on site i and spin a) and nfj = flafia 
(*fc<T = :^ Sfe e*''^*Si<j). The dispersion of the 
conduction band is e(k) and ef is the position 
of the localized level. The hybridization strength 
V is taken to be /c- independent, and finally U is 
the on-site Coulomb interaction strength between 
two /-electrons. Throughout this paper, the con- 
duction band will be described by a free (Bloch) 
density of states, po{E) = jf^^S{E — e(fc)), of 
semi-elliptic shape. Its width W = I sets the en- 
ergy scale, and its center of gravity the energy- 
zero: Tii = jf-^j:£{k) = 0. To obtain the single- 
electron Green's function for this model, we ap- 
ply dynamical mean-field theory (DMFT) in com- 
bination with the modified perturbation theory 
(MPT) to solve the associated impurity prob- 
lem. The DMFT approach becomes exact for 
Z — > oo (Z being the coordination number) and 
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represents a well-defined local approximation for 
finite dimensions (finite Z). The DMFT-MPT ap- 
proach recovers the high-energy features of the 
Green's function up to order (-Jr)- Furthermore, 
its low-energy behaviour is at least qualitatively 
correct. The investigation of ferromagnetism in 
the Hubbard model p2[ as well as the Mott- 
Hubbard transition ]2|, |23[| have shown that this 
method is able to predict qualitatively correct phe- 
nomena of strongly correlated systems. Numerical 
results are obtained using a standard DMFT algo- 
rithm which basically consists of a self-consistency 
loop as follows: Starting with an initial guess for 
the self-energy for the lattice, the conduction elec- 
tron bath of an associated impurity problem is de- 
fined using the self- consistency equation Then 
this impurity problem is solved by some means, 
and its self-energy extracted. The latter is then 
taken to be the lattice self-energy and a new impu- 
rity model is defined via the self-consistency equa- 
tion. This loop is iterated until self-consistency is 
achieved. More details on practical calculations 
can be found in [p| p4|. 

Next we need to specify how the impurity model 
is solved. For this we employ the modified per- 
turbation theory. This method is based on the 
following ansatz for the self-energy psl p6[: 
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(2) 



a^r and (3cr are introduced as parameters to be 
determined later. 'E'^'^'^\e) is the second-order 
contribution to perturbation theory around the 
Hartree-Fock solution ||2^. Equation (|[) can be 
understood as the simplest possible ansatz which 
can, on the one hand, reproduce the perturba- 
tional result in the limit [/ — > 0, and, on the other 
hand, recovers the atomic limit for appropriately 
chosen and Pa- p^ . 

Using the perturbation theory around the 
Hartree-Fock solution introduces an ambiguity 
into the calculation. Within the self-consistent 
Hartree-Fock calculation, one can either choose 
the chemical potential to be equivalent to the 
chemical potential of the full MPT calculation, or 
take it as parameter /i to be fitted to another phys- 
ically motivated constraint. In reference |^ the 
Luttinger theorem p8| , or equivalently the Friedel 
sum rule p9| , pO| , was used to determine fi. As 
discussed in Ref. we use the physically mo- 
tivated condition of identical impurity occupation 
numbers for the Hartree-Fock and the full calcu- 
lation (n'"/'^^^ — nl/^) to determine /i, which also 
allows for a consistent extension of the method to 
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Figure 1: Zero-temperature phase diagram for the 
periodic Anderson model without disorder from 
ijl^, J7 = 4 and V = 0.2. In the shaded region, the 
system is ferromagnetic, elsewhere paramagnetic. 



finite temperatures [pi] , |l?]] . Except for symmetric 
parameters this will lead to an approximate fulfill- 
ment of the Luttinger theorem only [^o| . Bearing 
in mind that for disordered systems, the Luttinger 
theorem does not apply, this should not be a deci- 
sive disadvantage for this study. A more detailed 
analysis of the different possibilities to determine 
fl is found in reference where the DMFT- 
MPT was applied to the single-band Hubbard 
model. Finally, the parameters and (3^ have 
to be determined. Instead of using the "atomic" 
limit of V — as was done for example in refer- 
ences ^ ^ , we make use of the moments of 
the spectral density. This procedure is described 
in detail in references |3^, ^ . The result not only 
fulfills the V = limit, but also recovers the high- 
energy behavior of the Green's function up to the 
order (-^). 

The results concerning ferromagnetism in the 
periodic Anderson model have been discussed else- 
where ^ ^ . For the following it is important 
to note that in the intermediate-valence regime, 
more precisely for —W/2 < Cf < —W/4, a ferro- 
magnetic solution with finite Tc exists for a range 
of electron densities. The T = phase diagram 
is plotted in Fig. |l|. In this area of the phase di- 
agram, the ferromagnetic solution shows typical 
features of a band-ferromagnet H, |j. The origin 
of the ferromagnetic order lies in the competition 
of kinetic and potential energy. Ferromagnetic or- 
der is stabilized by high values of the density of 
states (DOS) close to the lower band edge. The 
electron density needs to be chosen that the chem- 
ical potential lies in this region of large DOS. So 
although Stoner's theory does not capture the 
right physics, his criterion for the occurrence of 
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band-ferromagnetism {Up{fi) ^ 1) turns out to 
remain valid ^. 

Let us now turn to the problem of disorder. A 
general extension of model ([^) to include diagonal 
(on-site) disorder is 

H ->H + H^"^''^ (3) 

(4) 

In this model, the /-electron energy (et), the on- 
site hybridization (V) and the interaction strength 
([/) can deviate from the value denoted before by 
Ae/_i, AVi and AUi, respectively. The distribu- 
tion for each of these quantities can be defined 
by a probability distribution function P(Ae/_i), 
P{AV,) and P{AU^). 

A standard method to solve electron systems 
with disorder is the well-known coherent poten- 
tial approximation (CPA) |Tl|] . This method is 
considered the best single-site approximation and 
at least for one-particle properties, has proven to 
be remarkably successful It is known to be- 
come exact in a number of limiting cases, namely 
for small impurity concentration, small potential 
strengths, vanishing inter-site hopping, and as dis- 
cussed by ^ also for Z^^ . This last limiting 
case suggests that CPA results could be obtained 
by an alternative (DMFT-like) algorithm employ- 
ing a mapping onto an impurity model. This was 
shown to be the case and allows for a sys- 
tematic extension to include many-body interac- 
tions, which would be not possible within the stan- 
dard CPA procedure One has to bear 
in mind, however, that single-site approximations 
such as the one described here do have some limi- 
tations. A major limitation in regards to disorder 
is the inability to describe inhomogeneous, phase- 
separated systems which are discussed in the con- 
text of manganites jss). 

To calculate the single-electron Green function 
for a (diagonally) disordered system using DMFT, 
the algorithm described above for pure systems 
needs to be modified in the following way: In- 
stead of one, several impurity models need to be 
solved, one for each possible on-site configuration 
of the lattice model. The configurational aver- 
aging is then performed on the results for these 
impurity models, and the averaged self-energy ex- 
tracted and taken as lattice self-energy. This self- 
energy is fed into the self-consistency equation to 
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Figure 2: Curie temperatures for the PAM with 
disorder for [/ = 4, = 1.3, V = 0.2. The et 

and the disorder parameters are explained in the 
text. 



determine a new impurity bath function. This 
cycle is then iterated until self-consistency is ob- 
tained (i|, m. 

In the following, we want to present and dis- 
cuss numerical results for the periodic Anderson 
model with diagonal disorder. For simplicity, we 
have set all A[/j = and AVi = 0, including 
only disorder with respect to the /-level energy 
Cf. We further restrict ourselves to binary alloys, 
i.e. P(Aef,i) = p5{Aei) + (1 - p)(5(0). The in 
general rather complex definition of the disorder 
is reduced to two parameters: the difference of /- 
level positions of the two components Acf , and the 
concentration p. 

To examine the influence of disorder on ferro- 
magnetic order, we investigate two different sce- 
narios using this simple model of disorder: Select- 
ing U , V and the electron density so that 
a ferromagnetic solution is possible for a range of 
values of ef , we can choose ef and Aef so that one 
of the two alloy components would in a pure sys- 
tem be ferromagnetic, and the other not. We have 
taken ef = —0.4 and Acf = 0.2. These parame- 
ters for the two alloy components are indicated in 
the phase diagram (Fig. ^ as solid circles. The 
other scenario corresponds to alloying two differ- 
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Figure 3: Density of states A„{E) = 
-ilmGi*^(£') and Iml]^(£') for small inter- 
action strength U = 0.2, ef = -0.5, Aef = 0.2, 
V — 0.2 and n^*"*^ = 1.3 for the two pure limits 
p = and p = 1 as well as for p — 0.5 



ent ferromagnetic materials. It is realized by using 
Ef = —0.5 and Aet = 0.2. The resulting compo- 
nents are shown by the crosses in Fig. |l|. 

Let us start by looking at the first scenario: For 
p — the system is ferro-, for p ~ 1 paramagnetic. 
In Fig. 1^, we have plotted the Curie temperature 
as function of p. As one would expect, Tc decreases 
with increasing p, in the beginning linearly, and 
around p « 0.28, rather suddenly. With p > 0.29, 
the system is paramagnetic already for T — 0. 

The second scenario reveals more unexpected 
behaviour. For p = as well as for p = 1 (corre- 
sponding to the two crosses in Fig. |l|) the system is 
ferromagnetic and has finite Tc- In between these 
two pure limits, we find a strong reduction of the 
Curie temperature. Around p « 0.7, Tc even be- 
comes zero and the system is then paramagnetic. 
The reduction of Tc can be ascribed to the addi- 
tional quasiparticle damping induced by the dis- 
order. This effect of quasiparticle damping was 
already noticed in previous studies investigating 
the role of quasiparticle damping in the Hubbard 
model Q and in the PAM [|8| . The cited works 
did not involve disorder, but tested the influence 
of quasiparticle damping by comparing different 
approximation schemes some of which neglected 
quasiparticle damping completely. A similar re- 
duction of Tc due to disorder was also found for a 
ferromagnetic Kondo-lattice model with classical 
spins pq. 



To display the enhanced quasiparticle damping 
due to disorder, we have performed calculations 
for a weakly interacting (U = 0.2), paramagnetic 
system. In Fig. ||, the resulting f-density of states 
is plotted together with the imaginary part of the 
self-energy. These two quantities were calculated 
for p ~ Q, p ~ 0.5 and p = 1. Except for the 
smaller value of U, all other parameters were taken 
as in the lower panel of Fig. ||. For both pure sit- 
uations, the self-energy vanishes quadratically at 
the Fermi energy which is indicated by the respec- 
tive thin vertical line (note: the chemical poten- 
tial is shifted by varying p since we keep the total 
electron density n^'"*-' constant). In the alloyed 
compound, the imaginary part of the self-energy 
remains finite at the Fermi energy. And over a 
large energy range it is strongly enhanced com- 
pared to the pure limits. This strong disorder- 
induced quasiparticle damping is independent of 
the interaction strength, and can have a suppress- 
ing effect on ferromagnetism . 

Another effect is complementing the quasipar- 
ticle damping as mechanism to completely sup- 
press ferromagnetism as seen in Fig. |. As al- 
ready mentioned before, the occurrence of band- 
ferromagnetism is linked to a Stoner-like criterion 
requiring a large density of states at or close to the 
Fermi energy |j. For both pure cases {p — 
and p = 1), the Fermi energy lies within the charge 
excitation peak in the density of states as can be 
seen from Fig. ^ for the small U case. In the dis- 
ordered case, however, the structure of the DOS is 
dominated not by one, but two charge excitation 
peaks. The positions of these are given by the re- 
spective charge excitations in the pure limits. For 
intermediate values of p (p ~ 0.5), the Fermi en- 
ergy lies in between these two peaks, and the value 
of the density of states at this energy is relatively 
low. This should therefore reduce the tendency 
towards ferromagnetism. The special structure 
of the DOS for intermediate p leads therefore to 
further suppression of ferromagnetism. Whether 
(and for which p) Tc really vanishes, or not, now 
critically depends on the position of the Fermi en- 
ergy, and therefore on the electron density. 

To summarize, we have performed dynamical 
mean-field theory calculations for a periodic An- 
derson model (PAM) with disorder. The inclu- 
sion of disorder into these calculations follows a 
relatively simple and straightforward recipe p^ . 
Our numerical analysis shows how even 
small amounts of disorder can reduce the Curie 
temperature of a ferromagnetically ordered PAM 
in the intermediate-valence regime significantly. 
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Stronger disorder can lead to a complete suppres- 
sion of the ferromagnetic phase, even for a bi- 
nary alloy of two ferromagnetically ordered com- 
ponents. 
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